Controlled mechanical loading affects the osteocyte transcriptome in porcine trabecular bone in situ.

INTRODUCTION: Osteocytes modulate bone adaptation in response to mechanical stimuli imparted by the deforming bone tissue in which they are encased by communicating with osteoclasts and osteoblasts as well as other osteocytes in the lacuna-canalicular network through secreted cytokines and chemokines. Understanding the transcriptional response of osteocytes to mechanical stimulation in situ could identify new targets to inhibit bone loss or enhance bone formation in the presence of diseases like osteoporosis or metastatic cancer. We compared the mechanically regulated transcriptional response of osteocytes in trabecular bone following one or three days of controlled mechanical loading. METHODS: Porcine trabecular bone explants were cultured in a bioreactor for 48 h and subsequently loaded twice a day for one day or 3 days. RNA was isolated and sequenced, and the Tuxedo suite was used to identify differentially expressed genes and pathway analysis was conducted using Ingenuity Pathway Analysis (IPA). RESULTS: There were about 4000 differentially expressed genes following in situ culture relative to fresh bone. One hundred six genes were differentially expressed between the loaded and non-loaded groups following one day of loading compared to 913 genes after 3 d of loading. Only 45 of these were coincident between the two time points, indicating an evolving transcriptome. Clustering and principal component analysis indicated differences between the loaded and non-loaded groups after 3 d of loading. DISCUSSION: With sustained loading, there was a nine-fold increase in the number of differentially expressed genes, suggesting that osteocytes respond to loading through sequential activation of downstream genes in the same pathways. The differentially expressed genes were related to osteoarthritis, osteocyte, and chondrocyte signaling pathways. We noted that NFkB and TNF signaling are affected by early loading and this may drive downstream effects on the mechanobiological response. Moreover, these genes may regulate catabolic effects of mechanical disuse through their actions on pre-osteoclasts in the bone marrow niche.

Clinical outcomes of a prospective randomized comparison of bioreactance monitoring versus pulse-contour analysis in a stroke-volume based goal-directed fluid resuscitation protocol in brain-dead organ donors.

Eighty percent of brain-dead (BD) organ donors develop hypotension and are frequently hypovolemic. Fluid resuscitation in a BD donor is controversial. We have previously published our 4-h goal-directed stroke volume (SV)-based fluid resuscitation protocol which significantly decreased time on vasopressors and increased transplanting four or more organs. The SV was measured by pulse-contour analysis (PCA) or an esophageal doppler monitor, both of which are invasive. Thoracic bioreactance (BR) is a non-invasive portable technology that measures SV but has not been studied in BD donors. We performed a randomized prospective comparative study of BR versus PCA technology in our fluid resuscitation protocol in BD donors. Eighty-four donors (53.1%) were randomized to BR and 74 donors to PCA (46.8%). The two groups were well matched based on 24 demographic, social, and initial laboratory factors, without any significant differences between them. There was no difference in the intravenous fluid infused over the 4-h study period [BR 2271 ± 823 vs. PCA 2230 ± 962 mL; p = .77]. There was no difference in the time to wean off vasopressors [BR 108.8 ± 61.8 vs. PCA 150.0 ± 68 min p = .07], nor in the number of donors off vasopressors at the end of the protocol [BR 16 (28.6%) vs. PCA 15 (29.4%); p = .92]. There was no difference in the total number of organs transplanted per donor [BR 3.25 ± 1.77 vs. PCA 3.22 ± 1.75; p = .90], nor in any individual organ transplanted. BR was equivalent to PCA in clinical outcomes and provides a simple, non-invasive, portable technology to monitor fluid resuscitation in organ donors.

Quality control recommendations for RNASeq using FFPE samples based on pre-sequencing lab metrics and post-sequencing bioinformatics metrics.

BACKGROUND: Formalin-fixed, paraffin-embedded (FFPE) tissues have many advantages for identification of risk biomarkers, including wide availability and potential for extended follow-up endpoints. However, RNA derived from archival FFPE samples has limited quality. Here we identified parameters that determine which FFPE samples have the potential for successful RNA extraction, library preparation, and generation of usable RNAseq data. METHODS: We optimized library preparation protocols designed for use with FFPE samples using seven FFPE and Fresh Frozen replicate pairs, and tested optimized protocols using a study set of 130 FFPE biopsies from women with benign breast disease. Metrics from RNA extraction and preparation procedures were collected and compared with bioinformatics sequencing summary statistics. Finally, a decision tree model was built to learn the relationship between pre-sequencing lab metrics and qc pass/fail status as determined by bioinformatics metrics. RESULTS: Samples that failed bioinformatics qc tended to have low median sample-wise correlation within the cohort (Spearman correlation < 0.75), low number of reads mapped to gene regions (< 25 million), or low number of detectable genes (11,400 # of detected genes with TPM > 4). The median RNA concentration and pre-capture library Qubit values for qc failed samples were 18.9 ng/ul and 2.08 ng/ul respectively, which were significantly lower than those of qc pass samples (40.8 ng/ul and 5.82 ng/ul). We built a decision tree model based on input RNA concentration, input library qubit values, and achieved an F score of 0.848 in predicting QC status (pass/fail) of FFPE samples. CONCLUSIONS: We provide a bioinformatics quality control recommendation for FFPE samples from breast tissue by evaluating bioinformatic and sample metrics. Our results suggest a minimum concentration of 25 ng/ul FFPE-extracted RNA for library preparation and 1.7 ng/ul pre-capture library output to achieve adequate RNA-seq data for downstream bioinformatics analysis.

SARS-CoV-2 Variant Tracking and Mitigation During In-Person Learning at a Midwestern University in the 2020-2021 School Year.

Importance: The COVID-19 pandemic led many higher education institutions to close campuses during the 2020-2021 academic year. As campuses prepared for a return to in-person education, many institutions were mandating vaccines for students and considering the same for faculty and staff. Objective: To determine the association between vaccination coverage and the levels and spread of SARS-CoV-2, even in the presence of highly-transmissible variants and congregate living, at a midsized university in the US. Design, Setting, and Participants: This case series was conducted at a midsized Midwestern university during the spring 2021 semester. The university developed a saliva-based surveillance program capable of high-throughput SARS-CoV-2 polymerase chain reaction testing and genomic sequencing with the capacity to deliver results in less than 24 hours. On April 7, 2021, the university announced a vaccine requirement for all students for the fall 2021 semester and announced the same requirement for faculty and staff on May 20, 2021. The university hosted an onsite mass vaccination clinic using the 2-dose Pfizer-BioNTech vaccine during April 8 to 15 and April 29 to May 6, 2021. Data were analyzed for 14 894 individuals from the university population who were tested for COVID-19 on campus from January 6 to May 20, 2021. Main Outcomes and Measures: Positive SARS-CoV-2 diagnosis was confirmed by quantitative reverse transcription-polymerase chain reaction of saliva specimens, and variant identity was assessed by quantitative reverse transcription-polymerase chain reaction and next-generation sequencing of viral genomes. Results: Between January 6 and May 20, 2021, the university conducted 196 185 COVID-19 tests for 14 894 individuals and identified 1603 positive cases. Within those positive cases, 950 individuals (59.3%) were male, 644 (40.2%) were female, 1426 (89.0%) were students, and 1265 (78.9%) were aged 17 to 22 years. Among the 1603 positive cases, 687 were identified via polymerase chain reaction of saliva specimens. The Alpha (B.1.1.7) variant constituted 218 of the 446 total positives sequenced (48.9%). By May 20, 2021, 10 068 of 11 091 students (90.8%), 814 of 883 faculty (92.2%), and 2081 of 2890 staff (72.0%) were vaccinated. The 7-day rolling average of positive cases peaked at 37 cases on February 17 but declined to zero by May 14, 2021. The 7-day moving average of positive cases was inversely associated with cumulative vaccination coverage, with a statistically significant Pearson correlation coefficient of -0.57 (95% CI, -0.68 to -0.44). Conclusions and Relevance: This case series study elucidated the association of a robust vaccination program with a statistically significant decrease in positive COVID-19 cases among the study population even in the presence of highly transmissible variants and congregate living.

Performance of Three Tests for SARS-CoV-2 on a University Campus Estimated Jointly with Bayesian Latent Class Modeling.

Accurate tests for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) have been critical in efforts to control its spread. The accuracy of tests for SARS-CoV-2 has been assessed numerous times, usually in reference to a gold standard diagnosis. One major disadvantage of that approach is the possibility of error due to inaccuracy of the gold standard, which is especially problematic for evaluating testing in a real-world surveillance context. We used an alternative approach known as Bayesian latent class modeling (BLCM), which circumvents the need to designate a gold standard by simultaneously estimating the accuracy of multiple tests. We applied this technique to a collection of 1,716 tests of three types applied to 853 individuals on a university campus during a 1-week period in October 2020. We found that reverse transcriptase PCR (RT-PCR) testing of saliva samples performed at a campus facility had higher sensitivity (median, 92.3%; 95% credible interval [CrI], 73.2 to 99.6%) than RT-PCR testing of nasal samples performed at a commercial facility (median, 85.9%; 95% CrI, 54.7 to 99.4%). The reverse was true for specificity, although the specificity of saliva testing was still very high (median, 99.3%; 95% CrI, 98.3 to 99.9%). An antigen test was less sensitive and specific than both of the RT-PCR tests, although the sample sizes with this test were small and the statistical uncertainty was high. These results suggest that RT-PCR testing of saliva samples at a campus facility can be an effective basis for surveillance screening to prevent SARS-CoV-2 transmission in a university setting. IMPORTANCE Testing for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been vitally important during the COVID-19 pandemic. There are a variety of methods for testing for this virus, and it is important to understand their accuracy in choosing which one might be best suited for a given application. To estimate the accuracy of three different testing methods, we used a data set collected at a university that involved testing the same samples with multiple tests. Unlike most other estimates of test accuracy, we did not assume that one test was perfect but instead allowed for some degree of inaccuracy in all testing methods. We found that molecular tests performed on saliva samples at a university facility were similarly accurate as molecular tests performed on nasal samples at a commercial facility. An antigen test appeared somewhat less accurate than the molecular tests, but there was high uncertainty about that.

High-Throughput Genotyping of Common Chromosomal Inversions in the Afrotropical Malaria Mosquito Anopheles Funestus.

Polymorphic chromosomal inversions have been implicated in local adaptation. In anopheline mosquitoes, inversions also contribute to epidemiologically relevant phenotypes such as resting behavior. Progress in understanding these phenotypes and their mechanistic basis has been hindered because the only available method for inversion genotyping relies on traditional cytogenetic karyotyping, a rate-limiting and technically difficult approach that is possible only for the fraction of the adult female population at the correct gonotrophic stage. Here, we focus on an understudied malaria vector of major importance in sub-Saharan Africa, Anopheles funestus. We ascertain and validate tag single nucleotide polymorphisms (SNPs) using high throughput molecular assays that allow rapid inversion genotyping of the three most common An. funestus inversions at scale, overcoming the cytogenetic karyotyping barrier. These same inversions are the only available markers for distinguishing two An. funestus ecotypes that differ in indoor resting behavior, Folonzo and Kiribina. Our new inversion genotyping tools will facilitate studies of ecotypic differentiation in An. funestus and provide a means to improve our understanding of the roles of Folonzo and Kiribina in malaria transmission.

Inversion Genotyping in the Anopheles gambiae Complex Using High-Throughput Array and Sequencing Platforms.

Chromosomal inversion polymorphisms have special importance in the Anopheles gambiae complex of malaria vector mosquitoes, due to their role in local adaptation and range expansion. The study of inversions in natural populations is reliant on polytene chromosome analysis by expert cytogeneticists, a process that is limited by the rarity of trained specialists, low throughput, and restrictive sampling requirements. To overcome this barrier, we ascertained tag single nucleotide polymorphisms (SNPs) that are highly correlated with inversion status (inverted or standard orientation). We compared the performance of the tag SNPs using two alternative high throughput molecular genotyping approaches vs. traditional cytogenetic karyotyping of the same 960 individual An. gambiae and An. coluzzii mosquitoes sampled from Burkina Faso, West Africa. We show that both molecular approaches yield comparable results, and that either one performs as well or better than cytogenetics in terms of genotyping accuracy. Given the ability of molecular genotyping approaches to be conducted at scale and at relatively low cost without restriction on mosquito sex or developmental stage, molecular genotyping via tag SNPs has the potential to revitalize research into the role of chromosomal inversions in the behavior and ongoing adaptation of An. gambiae and An. coluzzii to environmental heterogeneities.

ColoType: a forty gene signature for consensus molecular subtyping of colorectal cancer tumors using whole-genome assay or targeted RNA-sequencing.

Colorectal cancer (CRC) tumors can be partitioned into four biologically distinct consensus molecular subtypes (CMS1-4) using gene expression. Evidence is accumulating that tumors in different subtypes are likely to respond differently to treatments. However, to date, there is no clinical diagnostic test for CMS subtyping. In this study, we used novel methodology in a multi-cohort training domain (n = 1,214) to develop the ColoType scores and classifier to predict CMS1-4 based on expression of 40 genes. In three validation cohorts (n = 1,744, in total) representing three distinct gene-expression measurement technologies, ColoType predicted gold-standard CMS subtypes with accuracies 0.90, 0.91, 0.88, respectively. To accommodate for potential intratumoral heterogeneity and tumors of mixed subtypes, ColoType was designed to report continuous scores measuring the prevalence of each of CMS1-4 in a tumor, in addition to specifying the most prevalent subtype. For analysis of clinical specimens, ColoType was also implemented with targeted RNA-sequencing (Illumina AmpliSeq). In a series of formalin-fixed, paraffin-embedded CRC samples (n = 49), ColoType by targeted RNA-sequencing agreed with subtypes predicted by two independent methods with accuracies 0.92, 0.82, respectively. With further validation, ColoType by targeted RNA-sequencing, may enable clinical application of CMS subtyping with widely-available and cost-effective technology.

Discovery and Structural Optimization of Acridones as Broad-Spectrum Antimalarials.

Malaria remains one of the deadliest diseases in the world today. Novel chemoprophylactic and chemotherapeutic antimalarials are needed to support the renewed eradication agenda. We have discovered a novel antimalarial acridone chemotype with dual-stage activity against both liver-stage and blood-stage malaria. Several lead compounds generated from structural optimization of a large library of novel acridones exhibit efficacy in the following systems: (1) picomolar inhibition of in vitro Plasmodium falciparum blood-stage growth against multidrug-resistant parasites; (2) curative efficacy after oral administration in an erythrocytic Plasmodium yoelii murine malaria model; (3) prevention of in vitro Plasmodium berghei sporozoite-induced development in human hepatocytes; and (4) protection of in vivo P. berghei sporozoite-induced infection in mice. This study offers the first account of liver-stage antimalarial activity in an acridone chemotype. Details of the design, chemistry, structure-activity relationships, safety, metabolic/pharmacokinetic studies, and mechanistic investigation are presented herein.

Draft Genome Sequences of Six Skin Isolates of Streptococcus pyogenes.

Whole-genome shotgun sequences and bottom-up assembly of contigs of six skin isolates of Streptococcus pyogenes, viz, NS88.3 (emm98.1), NS223 (emm91), NS455 (emm52), SS1448 (emm86.2), SS1572 (emm223), and SS1574 (emm224), are presented here. All contigs were annotated, and the gene arrangements and the inferred proteins were consistent with a pattern D classification.

High-quality genetic mapping with ddRADseq in the non-model tree Quercus rubra.

BACKGROUND: Restriction site associated DNA sequencing (RADseq) has the potential to be a broadly applicable, low-cost approach for high-quality genetic linkage mapping in forest trees lacking a reference genome. The statistical inference of linear order must be as accurate as possible for the correct ordering of sequence scaffolds and contigs to chromosomal locations. Accurate maps also facilitate the discovery of chromosome segments containing allelic variants conferring resistance to the biotic and abiotic stresses that threaten forest trees worldwide. We used ddRADseq for genetic mapping in the tree Quercus rubra, with an approach optimized to produce a high-quality map. Our study design also enabled us to model the results we would have obtained with less depth of coverage. RESULTS: Our sequencing design produced a high sequencing depth in the parents (248×) and a moderate sequencing depth (15×) in the progeny. The digital normalization method of generating a de novo reference and the SAMtools SNP variant caller yielded the most SNP calls (78,725). The major drivers of map inflation were multiple SNPs located within the same sequence (77% of SNPs called). The highest quality map was generated with a low level of missing data (5%) and a genome-wide threshold of 0.025 for deviation from Mendelian expectation. The final map included 849 SNP markers (1.8% of the 78,725 SNPs called). Downsampling the individual FASTQ files to model lower depth of coverage revealed that sequencing the progeny using 96 samples per lane would have yielded too few SNP markers to generate a map, even if we had sequenced the parents at depth 248×. CONCLUSIONS: The ddRADseq technology produced enough high-quality SNP markers to make a moderately dense, high-quality map. The success of this project was due to high depth of coverage of the parents, moderate depth of coverage of the progeny, a good framework map, an optimized bioinformatics pipeline, and rigorous premapping filters. The ddRADseq approach is useful for the construction of high-quality genetic maps in organisms lacking a reference genome if the parents and progeny are sequenced at sufficient depth. Technical improvements in reduced representation sequencing (RRS) approaches are needed to reduce the amount of missing data.

Genome-wide profiling of 24 hr diel rhythmicity in the water flea, Daphnia pulex: network analysis reveals rhythmic gene expression and enhances functional gene annotation.

BACKGROUND: Marine and freshwater zooplankton exhibit daily rhythmic patterns of behavior and physiology which may be regulated directly by the light:dark (LD) cycle and/or a molecular circadian clock. One of the best-studied zooplankton taxa, the freshwater crustacean Daphnia, has a 24 h diel vertical migration (DVM) behavior whereby the organism travels up and down through the water column daily. DVM plays a critical role in resource tracking and the behavioral avoidance of predators and damaging ultraviolet radiation. However, there is little information at the transcriptional level linking the expression patterns of genes to the rhythmic physiology/behavior of Daphnia. RESULTS: Here we analyzed genome-wide temporal transcriptional patterns from Daphnia pulex collected over a 44 h time period under a 12:12 LD cycle (diel) conditions using a cosine-fitting algorithm. We used a comprehensive network modeling and analysis approach to identify novel co-regulated rhythmic genes that have similar network topological properties and functional annotations as rhythmic genes identified by the cosine-fitting analyses. Furthermore, we used the network approach to predict with high accuracy novel gene-function associations, thus enhancing current functional annotations available for genes in this ecologically relevant model species. Our results reveal that genes in many functional groupings exhibit 24 h rhythms in their expression patterns under diel conditions. We highlight the rhythmic expression of immunity, oxidative detoxification, and sensory process genes. We discuss differences in the chronobiology of D. pulex from other well-characterized terrestrial arthropods. CONCLUSIONS: This research adds to a growing body of literature suggesting the genetic mechanisms governing rhythmicity in crustaceans may be divergent from other arthropod lineages including insects. Lastly, these results highlight the power of using a network analysis approach to identify differential gene expression and provide novel functional annotation.

An SNP-based second-generation genetic map of Daphnia magna and its application to QTL analysis of phenotypic traits.

BACKGROUND: Although Daphnia is increasingly recognized as a model for ecological genomics and biomedical research, there is, as of yet, no high-resolution genetic map for the genus. Such a map would provide an important tool for mapping phenotypes and assembling the genome. Here we estimate the genome size of Daphnia magna and describe the construction of an SNP array based linkage map. We then test the suitability of the map for life history and behavioural trait mapping. The two parent genotypes used to produce the map derived from D. magna populations with and without fish predation, respectively and are therefore expected to show divergent behaviour and life-histories. RESULTS: Using flow cytometry we estimated the genome size of D. magna to be about 238 mb. We developed an SNP array tailored to type SNPs in a D. magna F2 panel and used it to construct a D. magna linkage map, which included 1,324 informative markers. The map produced ten linkage groups ranging from 108.9 to 203.6 cM, with an average distance between markers of 1.13 cM and a total map length of 1,483.6 cM (Kosambi corrected). The physical length per cM is estimated to be 160 kb. Mapping infertility genes, life history traits and behavioural traits on this map revealed several significant QTL peaks and showed a complex pattern of underlying genetics, with different traits showing strongly different genetic architectures. CONCLUSIONS: The new linkage map of D. magna constructed here allowed us to characterize genetic differences among parent genotypes from populations with ecological differences. The QTL effect plots are partially consistent with our expectation of local adaptation under contrasting predation regimes. Furthermore, the new genetic map will be an important tool for the Daphnia research community and will contribute to the physical map of the D. magna genome project and the further mapping of phenotypic traits. The clones used to produce the linkage map are maintained in a stock collection and can be used for mapping QTLs of traits that show variance among the F2 clones.

Genome-wide polymorphisms and development of a microarray platform to detect genetic variations in Plasmodium yoelii.

The rodent malaria parasite Plasmodium yoelii is an important model for studying malaria immunity and pathogenesis. One approach for studying malaria disease phenotypes is genetic mapping, which requires typing a large number of genetic markers from multiple parasite strains and/or progeny from genetic crosses. Hundreds of microsatellite (MS) markers have been developed to genotype the P. yoelii genome; however, typing a large number of MS markers can be labor intensive, time consuming, and expensive. Thus, development of high-throughput genotyping tools such as DNA microarrays that enable rapid and accurate large-scale genotyping of the malaria parasite will be highly desirable. In this study, we sequenced the genomes of two P. yoelii strains (33X and N67) and obtained a large number of single nucleotide polymorphisms (SNPs). Based on the SNPs obtained, we designed sets of oligonucleotide probes to develop a microarray that could interrogate ∼11,000 SNPs across the 14 chromosomes of the parasite in a single hybridization. Results from hybridizations of DNA samples of five P. yoelii strains or cloned lines (17XNL, YM, 33X, N67 and N67C) and two progeny from a genetic cross (N67×17XNL) to the microarray showed that the array had a high call rate (∼97%) and accuracy (99.9%) in calling SNPs, providing a simple and reliable tool for typing the P. yoelii genome. Our data show that the P. yoelii genome is highly polymorphic, although isogenic pairs of parasites were also detected. Additionally, our results indicate that the 33X parasite is a progeny of 17XNL (or YM) and an unknown parasite. The highly accurate and reliable microarray developed in this study will greatly facilitate our ability to study the genetic basis of important traits and the disease it causes.

Naloxone-precipitated conditioned taste aversions in morphine-dependent Fischer (F344) and Lewis rat strains.

The Fischer 344 (F344) and Lewis (LEW) rat strains are genetically divergent populations that are used to study the effects of and responses to drugs of abuse. In this context, LEW rats display faster acquisition of drug self-administration than F344 rats. Interestingly, these strains have also been reported to differ in their somatic responses to morphine withdrawal. To address possible strain differences in the affective response to withdrawal, the present study assessed the ability of naloxone-precipitated withdrawal from morphine to induce conditioned taste aversions in male F344 and LEW rats. Specifically, subjects from each of these strains were given chronic morphine to induce dependence and then given access to a novel saccharin solution followed by naloxone. These pairings were given every fourth day for a total of two conditioning trials after which subjects were given access to saccharin but without naloxone administration to assess extinction of the naloxone-induced aversion. Behavioral assays of withdrawal were also performed after each naloxone administration. Both F344 and LEW subjects acquired aversions to the naloxone-associated taste with no significant differences in the rate of acquisition of the aversions. Differences did appear during extinction with LEW animals extinguishing the taste aversion significantly faster than F344 animals. The data were discussed in terms of the relative strength of the affective responses during withdrawal and the role of such responses to drug use and abuse.